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Fig. 3.6 Mean free path of radiation in NO»; radiation the same temperature 
as the NOp. X= 1/K! P(NOp) 
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Fig. 3.7 Mean free path of radiation in NO»; the radiation temperature above 
10,000°K 


~ 85 - 


{e) 
(o) 
epeseisee rH ro 
rH N 
tt + i 
ane root H Ha : 4 ce 
it siieerrti Searant cet Ht ania) 
: [e) 
pei Ferree i HER ! | 
~~ 
2 3 
2) + 
= tl 
Py 
fo) Se 
: R i 
= ras) 
So 
= on 
4 Sy 
HG Sx 
‘ Or — 
H ; ~ & on 1 
a se ! 
‘iy a | 2 
+ ro! ue 
it } uw ° | 
Et ow 
EH a I 
: = 
4 fo) amd 
t g 2. 
t+ ° 
7 2 
HH <q 
4 8 
Mm ae 
H me inp) 
: % 
H 8 Fy 
BRepeaaGuate, N 
a 9 
zat 2) 


oO oO °o 
ad wo ve) Tr 


(WO 71) Mtv 4O LN319143909 NOLLaYOSaY 


6) 


9 


9 


OABM JSBIG OY} Jo AT0}STY UIT} [[B-10AQ OTE “STA 


6) 


6) 


3 


(SYALSW) SNIGVY YOOHS 
92 og Zl 


(SGNOOJS) AWIL 
el 


= 96 = 


below 5000°K. Of course the air has no NOs» at all before it enters the 
shock and so some time is required to form NO, after the air is heated by 
the shock. Formation of NO, in heated air under laboratory conditions has 
been studied by Daniels. Presumably the mechanism for formation is in 
two steps 


(a) N, + 0, + 2 NO 


(b) 2NO +0, + 2 NO, 


Reaction (b)!5 is extremely fast at all temperatures above 600°K and will 
always maintain equilibrium between NO, Og and NO»g. The rate of reaction 
(a) is strongly temperature dependent and so effectively one can say that it 
takes place above a certain temperature and does not occur below this tem- 
perature. We have selected 2000°K for this temperature from data furnished 
by Daniels.!¢ 

The NO» zone is therefore limited to the air which was initially within 
100 meters, since the shock temperature is 2000°K at this distance. 

Last Stage of the Cooling. At the time the core starts radiating most 
strongly, the shock front is most opaque and therefore the radiated energy 


15. Hans J. Schumacker, Chemische Gasreactionen, T. Steinkopff, Dresden 
und Leipzig, 1938. 


16. Farrington Daniels, private communication. 
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Air at the shock front when the temperature is 2000°K has a pressure 
of about 10 atm Oo and 40 atm No. 
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dt 


The rate of formation of NO is 0.019 atm/msec while the equilibrium 
partial pressure is 0.4 atm. The shock front moves about a meter per 
millisecond and so obviously the NO, equilibrium does not have time to 
be established. A similar calculation shows that at 2500°K the equilib- 
rium is easily established within a millisecond. The indicated NO» 
"cutoff" is thus below 2500°K and we have arbitrarily taken 2000°K as 
the limit. 


8.5 -22.640 
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10 x 400 = 400 x 107 1°33 = 18.7 atm/sec 
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is essentially all absorbed behind the front. There are about two optical 
depths for high temperature radiation due to the NO, alone, and at this time 
Oy in its lowest state is also very effective since the radiation has a tem- 
perature of about 40,000°K. We have not followed the decrease in opacity 

of the absorption layer as it expands, but by the time the shock front 

reaches 300 meters there is certainly negligible absorption of 20,000°K radi- 
ation by all outside layers. The ball of fire finally radiates as a hot mass 
of gas without self absorption. At the very last, only NO» can radiate visible 
light and so the ball of fire darkens as it drops to 2000°. 


3.6 The Ball of Fire 


The last two sections have described the effects of radiation in the 
formation and growth of the air blast. Except for the noted omissions,!” 
this description is qualitatively complete for air blasts of this size. Other 
phenomena of interest have to do with the transmission of radiation energy 
ahead of the shock. Insofar as the description of the last two sections gives 
a complete model of the air blast, we can use it to determine the amount of 
energy transmitted through the shock front as a function of time. We also 
get a description of the quality (e.g., temperature) of the radiation, and the 
size and appearance of the ball of fire. 

In the early phase, as long as the front is itself radiating strongly, 
there is a well-defined, spherical ball of fire. During this time the inner 
radiation front is, of course, not to be seen. When the shock front can no 
longer radiate, the hot core inside becomes more visible as the molecular 
absorbing layer becomes more tenuous. In the final stage, at the time of 
the second maximum and later, there is a hot core radiating, surrounded by 
a cooler fringe which is also radiating. At this time the ball has little ab- 
sorption for its own radiation and is no longer a well-defined sphere with a 
sharp boundary. 

The amount of radiant energy escaping per unit time into the air 
ahead of the shock front as a function of time is given in Fig. 3.10. This 
figure is largely schematic and has been obtained by methods which will be 
discussed below. 6 

The unit © (used in Fig. 3.10) is 1.35 x 10 ergs/cm2/sec, the flux of 
radiant energy incident into the earth's outer atmosphere from the sun. D 
is measured in meters. 


17. Effects of gamma radiation, bomb materials, turbulence. 


- 98 - 


sniper yooys sa oinyeredure}z yooys I1°S ‘sla 


SNIdVY 


i 


TTT 


(Ho) t 


- 102 - 


Chapter 6 


THE IBM SOLUTION OF THE BLAST WAVE PROBLEM 
by K. Fuchs 


6.1 Introduction 


The discussion in the preceding chapters shows that the problem of 
the propagation of the blast wave from a nuclear explosion is quite compli- 
cated. Even if we disregard any transport of radiation, appreciable compli- 
cations arise in view of the large variations in temperature and entropy. 
We may roughly divide the range of shock temperatures into two regions. 
First, the region from about one million degrees to about three thousand 
degrees absolute; here dissociation of molecules and ionization of the atoms 
takes place; consequently y — 1 is fairly small but varies with temperature. 
Second, below 3000°; here y is less variable and approaches eventually the 
value 1.4 in normal conditions. 

The temperature varies also, of course, along an adiabatic. However, 
the total variation between the shock pressure and 1 atm is not excessive 
— about a factor 2 at shock temperature of 3000°K and slightly more than a 
factor 10 for a shock temperature of 1,000,000°. The effect of decreasing 
temperature along the adiabatic on the degree of ionization or dissociation 
is partly balanced by the decrease in density. For this reason the variation 
of y along an adiabatic is not as pronounced as one might expect from the 
change in temperature. Therefore qualitative statements made about the 
conditions at the shock front hold to a large degree also for the subsequent 
expansion behind the shock. 

_A temperature of 3000°K is reached in the shock when the shock pres- 
sure is about 80 atm. (For an energy release of 10,000 tons of TNT the 
shock radius is then 80 meters.) We are, however, more interested in the 
pressure region from about 1 atm down (corresponding to shock radii of 
500 meters or more). It was felt that the exact energy distribution at this 
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early stage could not have an appreciable effect on the pressure distribution - 
in the later stages as long as the energy distribution is at least roughly 
correct. It was decided to start an IBM calculation at this point and to 
calculate the initial conditions for this instant by means of the approxima- 
tions developed in the preceding chapters. The approximations made in- 
clude (for details see Section 6.2): 

(1) An approximate treatment of the isothermal sphere. In actual 
fact the isothermal sphere at this late stage has no significant influence on 
the propagation of the shock. However, in the first instance we are inter- 
ested in the isothermal sphere as such; in the second instance the iso- 
thermal sphere greatly facilitates numerical computations, since it reduces 
the total range of temperatures and entropies which exist at any given mo- 
ment, and eliminates the singularity at the centre. 

(2) y was assumed to be constant and an average value 1.25 was as- 
sumed. This is probably the least satisfactory of the assumptions made, 
but in view of the fact that we did not require a very accurate estimate of 
the initial conditions, and that a calculation with variable y is hardly feasible 
without a great amount of computation, this assumption seemed justified. 

(3) y—1 was assumed small. This assumption is not essential, but 
the error introduced thereby is small and it has the advantage that the 
isothermal sphere can be included as an integral part of the calculation. 

For the IBM run it is of great advantage if the variation of the pres- 
sure along an adiabatic is a simple function of the density, though the 
variation from one adiabatic to another may be given in numerical form. 

In this case the adiabatic of each mass point is given by one or two con- 5 
stants and we require only a table of these constants as functions of one 
variable. 

It would be difficult to find a simple presentation of the equation of 
state covering the whole region of adiabatics between the shock pressure 
curve and normal pressure. However, with the limitation to shock pressure 
below 80 atm the demand of a simple representation becomes feasible. In 
particular this is true for the adiabatics which start at shock pressures 
below 80 atm, since then no ionization or dissociation occurs. We require 
also, of course, the adiabatics of the inner mass points, which were shocked 
by stronger shocks. However we require only the tail end of these adia- 
batics. Furthermore, the highest entropies are eliminated by the equaliza- 
tion of entropy inside the isothermal sphere. The equation of state which 
has been used is discussed in Chapter 7. The requirement of a simple _ 
equation of state is the principal reason for starting the IBM run at such a 
comparatively late stage. 

Although radiation transport has been taken into account insofar as it 
is responsible for the formation of the isothermal sphere, no allowance has 
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been made for the radiation transport from the isothermal sphere into the 
region in which NO» is formed. As explained in Chapter 3, this transport 
of energy’ becomes important when the shock radius has reached about 100 
meters and it should affect the propagation of the shock shortly thereafter. 
The opacity data required for the purpose of calculating this transport are 
not sufficiently well known.! In neglecting the radiation transport altogether 
we are pessimistic, since it is of advantage to have the energy close to the 
shock front. Then the shock pressure decreases less rapidly than it would 
otherwise. The increased shock pressure would naturally lead to a greater 
degree of dissipation of energy by the shock, so that at larger distances the 
shock pressure might drop again more rapidly, and at sufficiently large 
distances the effect of the radiation transport on the shock pressure would 
be reversed. At present we are not in a position to make any definite 
statement about this possibility. 

The IBM run was intended for an energy release of 10,000 tons of 
TNT. Owing to some unfortunate circumstances related in Section 6.3, no 
definite energy can be attributed to the run throughout its whole history. 

At sufficiently large distances the energy should be assumed to be 13, 000 
tons. Other energies can, of course, be obtained by the usual scaling laws. 


6.2 The Initial Conditions of the IBM Run 


The initial conditions of the IBM run were prepared by Hirschfelder 
and Magee. The principal data are summarized below without going into the 
details of the calculation. All data are for an energy release of 10, 000 
tons of TNT. 


6.2.1 The Isothermal Sphere? 


The Lagrangian coordinate of the radiation front is given in terms of 
the shock radius by Eq. 3.64 of Chapter 3. It can be rewritten in the form 


0.7326 1/3 


ry = (11.85 ¥ ~ 601.5 100 cm (6.1) 


Here both rg and Y are given in centimeters. It was convenient to choose 


1. Chapter 3 contains a discussion of opacity data which postdates this 
statement. 


2. Chapter 3 values are revised figures. Notations in equations used here 
are in the main those of Bethe in Chapter 4. 
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a simple value of Y/ tr and the value 4 was chosen, which corresponds to 
a shock pressure very near to 80 atm. Then 


7) = 1997 cm, Y = 7987 cm, Y¥/t5 = 4 (6.2) 


The actual radius Ro of the isothermal sphere is obtained from the con- 
servation of mass. Since we assume constant density p and constant pres- 
sure in the isothermal sphere one has 


(6.3) 


The value rg varies very slowly after the shock radius has reached about 
10 meters, and the effect of the isothermal sphere on the shock is negli- 
gible a short while thereafter. Beyond a shock radius of 80 meters rg 
varies very slowly indeed. 

It is a good assumption to assume that Yg is constant. The actual 
radius Rp then varies in accordance with Eq. 6.3 only because the air in 
the isothermal sphere expands. The initial value of Ro is 60 meters. 

The initial condition of the isothermal sphere is given by the following 
quantities: 


Temperature = 49, 000°K 

Pressure = 37.0 atm 

Density = 0.0392 x normal density 
Entropy As/R = 85 


Internal energy and enthalpy E/R = 1.487 x 10°, H/R = 1.782 x 10° 
Ro = 60.23 meters 


ro = 19.97 meters 


These data were obtained from a calculation indicated below. 


6.2.2 Initial Pressure and Density Distribution 


It has been shown by Bethe in Chapter 4, Section 5, that the isothermal 
sphere can be treated on the assumption of small y —1. The small y-l1 
approximation has been checked for a point source solution (see Chapter 4, 
Sections 3 and 4) and it was found satisfactory in the region in which we 
are interested. Since the solution which we require is in any case very 
close to the point source solution except in the neighborhood of the isothermal 
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sphere, the error of the small (y — 1) approximation is small. 

We shall not go into the details of the calculation. The analysis is 
rather involved, but the lines along which it proceeds are sufficiently indi- 
cated in Bethe's Section 4.5. The resulting equations had previously been 
evaluated for two values of y in order to see how sensitive they were to a 
change in y. The values chosen were y = 1.2 and y = 1.3. The values for 
y = 1.25 were then obtained by interpolation. In this way we find the initial 
pressure and density distribution as well as the velocities and the Eulerian 
coordinates of all mass points. 

Some adjustments had to be made on the intial conditions. A minor 
adjustment arose from the fact that the small y — 1 treatment of the iso- 
thermal sphere does not agree exactly with the treatment of Hirschfelder and 
Magee. The error, which may be due to either method, may be seen from 
the values of Ro/Y; the small y = 1 treatment gives Ro/Y = 0.850 compared 
to 0.761 by means of the other method. The latter value was assumed to be 
more reliable. 

Furthermore, at the start of the IBM calculation the value of y at the 
shock front is larger than 1.25; instead of a compression ratio in the shock 
Ps/Po = 9, as would be expected for y = 1.25, the value obtained from the 
correct Hugoniot curve for a pressure of 77.25 atm is Pg/Po = 7.24. The 
density contour was, therefore, adjusted to give the correct compression 
ratio at the shock, and the correct radius of the isothermal sphere. This 
required also an adjustment in the Eulerian coordinates R, since it is essen- 
tial that the initial conditions satisfy the equation of continuity: 


— (6.4) 


where r is the Lagrangian coordinate. 
The initial velocities were then calculated directly from the equation 


( ne 
. Lt +\— 
Y Y 


Qe re 
+ yp ears 
1+y ¥ 


(6.5) 


which follows from the small y — 1 approximation. Here the correct shock 
velocity Y for the given shock pressure was used. 
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6.3 The Total Energy : 


The methods used to establish suitable initial conditions are in parts 
somewhat arbitrary. For this reason it is not surprising that the total 
energy corresponding to these conditions turned out to deviate appreciably 
from the assumed value of 10,000 tons of TNT. Unfortunately, the energy 
was recalculated from the initial conditions only after the IBM run had been 
completed. It was then found that the total energy was 13,500 tons TNT. 

Since the initial shock pressure of 77.25 atm at the initial shock radius 
of 79.9 meters corresponds to an energy release of 10,000 tons, we have no 
remedy for the discrepancy. All that can be said is that the shock pressure 
vs distance curve corresponds to 10,000 tons up to 80 meter shock radius 
and to 13,500 tons for large radii. For intermediate radii it should slowly 
change between these values. 

Actually the discrepancy is slightly less. A check of the total energy 
at a shock radius of 2000 meters gave only 13,100 tons. The "loss" of 
400 tons is entirely due to errors of the IBM run and is of the order of 
magnitude to be expected from this source. 

For most purposes the total energy in the IBM run should be assumed 
to be about 13, 000 tons, except at small shock radii, where 10,000 tons is 
more appropriate. 


6.4 The IBM Run 


The hydrodynamical equations are 


2 2 
8 RR ap 
Py zt = 0 (6.6) 
One | ge oF 
2 
Po oR” aR 
grr oa eo 
r 


In addition we have the equation of the adiabatics which were put in the 
form 


: » \8 
Cae (] + nf (6.8) 
~ 182 - : 


where g and 7 are numerically known functions of the entropy. 

The boundary conditions are that the velocity dR/dt vanish at the cen- 
ter and that at the shock radius the Hugoniot conditions be satisfied. The 
latter determine also the entropy of any mass point as it passes through 
the shock front. The entropy is assumed to remain constant, so that es- 
sentially g and 7 are given functions of the Lagrange variable r. 

The method employed first in solving the system of partial differential 
equations by means of the IBM machines was suitable as long as the shock 
pressure differed appreciably from 1 atm, but it became erratic as the 
overpressure became small. The method was therefore changed, so as to 
calculate changes in density and pressure rather than their absolute values. 
This change of procedure quickly suppressed the erratic behavior of the 
pressure. 

The run was continued until the shock radius had reached a value of 
6.270 meters. At that instant the overpressure in the shock was 0.0251 
atm. The positive pulse was 290 meters long and the negative pulse 760 
meters. Since the further propagation of the shock is influenced at these 
low overpressures only by the positive pulse, the approximations on which. 
the semi-acoustic theory of Chapter 5 are based are well satisfied. They 
are (1) that the overpressure be small compared to 1 atm, and (2) that the 
length of the pressure pulse be small compared to the shock radius. Even 
the application of the semi-acoustic theory to the negative phase is not bad. 
Hence, the IBM run was discontinued and the semi-acoustic theory was used 
for the purpose of continuation. 


6.5 Results 


The shock pressure as a function of the distance of the shock front 
from the center of the explosion is shown in Fig. 6.1. 

In this graph all data have been collected from the previous chapters, 
as well as from calculations which do not appear in this volume. From a 
shock radius of 10 to 80 meters the similarity solution has been used. The 
dotted lines show upper and lower limits for the effect of the bomb material 
on the propagation of the shock, as calculated for the Trinity test (sealed 
down to 10,000 tons of TNT). In this region the curve corresponds to an 
energy release of 10,000 tons. 

For shock radii from 80 to 6300 meters, the shock pressures are ob- 
tained from the IBM run. Here the total energy is between 10,000 and 
13,000 tons, the upper value being correct at sufficiently large distances. 
The time of arrival of the shock is indicated at various points. T = 0 is 
the start of the IBM run which was 0.012 sec after the explosion. 
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Beyond a radius of 6300 meters up to 67,000 meters, the semi-acoustic 
theory of Chapter 5 has been used. 

A number of curves showing the pressure at a fixed distance as func- 
tion of time are shown in Figs. 6.2 to 6.4. Graphs of the pressure vs 
distance at a fixed time are shown in Figs. 6.5 to 6.13. 

The duration of the positive phase of the pulse as function of the shock 
pressure is given in Fig. 6.14. 

Finally, there are shown in Fig. 6.15 the positive impulse I+ and the 
fraction of the total energy which is left in the blast as functions of the 
shock pressure divided by the normal pressure. The latter is independent 
of the energy release. The positive impulse has been scaled to an energy 
release of 40,000 tons in free air (or 20,000 tons on the ground) for the 
purpose of comparison with observations at Trinity. 


6.6 Comparison with TNT Explosion. Efficiency of Nuclear Explosion 
TS EE i chenicy OF Nuclear Explosion 


One purpose of the IBM run was to find the efficiency of a nuclear 
explosion compared to an explosion from an equivalent charge of TNT. In 
the nuclear explosion a greater amount of energy is used for the purpose 
of heating the air near the center of the explosion to high temperatures. 

A large fraction of this energy is useless for the propagation of the shock. 

Since no comparable IBM run exists for a TNT explosion, we compared 
the results for the nuclear explosion with experimental data. For this 
purpose, the experimental curve prepared by Hirschfelder, Littler and 
Sheard* was used. It is based on experimental data for charges fired on 
the ground, and for shock pressures in the range from 15 to 2 pounds per 
square inch. The charges varied from 67 to 550 pounds. The curve is 
given by the expression 


ADE Eee gee (6.9) 


Pe a (6.10) 


Here Y is the shock radius in feet; w the weight of the charge in pounds; 
and Ap the overpressure in psi. 


4. Hirschfelder, Littler and Sheard, Estimated Blast Pressures from TNT 
charges of 2 to 10,000 tons, Los Alamos Scientific Laboratory Report 
LA~316, June 25, 1945 (classified). 
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